Complex angular momentum in black hole physics and quasinormal modes 
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By using the complex angular momentum approach, we prove that the quasinormal mode complex 
frequencies of the Schwarzschild black hole are Breit-Wigner type resonances generated by a family 
of surface waves propagating close to the unstable circular photon (graviton) orbit at r = 3A/. 
Furthermore, because each surface wave is associated with a given Regge pole of the S-matrix, we 
can construct semiclassically the spectrum of the quasinormal-mode complex frequencies from Regge 
trajectories. The notion of surface wave orbiting around black holes thus appears as a fundamental 
concept which could be profitably introduced in various areas of black hole physics in connection 
with the complex angular momentum approach. 
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I. INTRODUCTION 

Since the pioneering work of Watson dealing with 
the propagation and diffraction of radio waves around the 
Earth, the complex angular momentum (CAM) method 
has been extensively used in several domains of scatter- 
ing theory (see the monographs of Newton [2| and of 
Nussenzveig 3j and references therein for various applica- 
tions in quantum mechanics, nuclear physics, electromag- 
netism, optics, acoustics and seismology). The success of 
the CAM method is due to its ability to provide a clear 
description of a given scattering problem by extracting 
the physical information (linked to the geometrical and 
diffractive aspects of the scattering process) which is hid- 
den in partial-wave representations. 

The CAM method was first used in gravitational 
physics by Chandrasekar and Ferrari £| in their study of 
nonradial oscillations of relativistic stars. They used the 
theory of Regge poles Q to determine the flow of grav- 
itational energy through the star. The general frame- 
work for the CAM description of Schwarzschild black 
hole scattering was developed by Andersson and Thylwe 
, which Andersson |(J then used to interpret the black 
hole glory. An important concept, which is naturally 
present in the CAM point of view, is that of a sur- 
face wave, which allows a description of the diffractive 
effects of scattering. Andersson established that the 
surface waves orbiting around the Schwarzschild black 
hole of mass M propagate close to the unstable pho- 
ton orbit at r — 3M. Yet except for these articles, 
the CAM/surf ace- wave approach to resonant scattering 



in black hole physics has been neglected in favor of the 
quasi-normal modes (QNMs), in which the dynamical re- 
sponse to an external perturbation is explained in terms 
of resonant frequencies. For recent reviews of the QNM 
approach see Refs. and0 as well as chapter 4 of Ref. 0. 
An introduction to black hole scattering can be found in 

In this paper we will establish the connection be- 
tween Andersson's surface waves and the QNMs. More 
precisely, we shall prove here that the QNM complex 
frequencies of the Schwarzschild black hole are Breit- 
Wigner-type resonances generated by the surface waves. 
Moreover, because each surface wave is associated with 
a given Regge pole of the S'-matrix, we can construct 
the spectrum of the QNM complex frequencies from the 
Regge trajectories, i.e., from the curves traced out in the 
CAM plane by the Regge poles as a function of the fre- 
quency. 

As early as 1972, it was suggested by Goebel that 
the black hole normal modes could be interpreted in 
terms of gravitational waves in spiral orbits close to the 
unstable photon orbit at r = 3M, which decay by radi- 
ating away energy. In the present paper, using the CAM 
approach, we establish this appealing and physically in- 
tuitive picture on a rigorous basis. To conclude, we also 
provide a framework for future developments. 



II. FROM CAM TO QNM 
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We first consider the scattering of a monochromatic 
scalar wave, with time dependence exp(— iuot), by the 
Schwarzschild black hole of mass M. The correspond- 
ing scattering amplitude can be written (see, e.g. 0) 
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f(cj,6) = — V(2£ + 1) (S e (uj) - 1) P e (cos9) 
liui L — ' 



(1) 



i=0 



where £ is the ordinary angular momentum index, Pg are 
the usual Legendre polynomials and Si are the diagonal 
elements of the ^-matrix. For a given angular momentum 
£, the coefficient St is obtained from the partial wave 
solution $f of the following problem: 
(i) <&i satisfies the Schrodinger-type equation 



dri 



[u, 2 - V(r)] $ t = 0. 



(2) 



Eq. J2J) is obtained from the scalar wave equation af- 
ter separating variables and is called the Regge- Wheeler 
equation. Here r is the standard radial Schwarzschild 
coordinate and r* = r + 2Mln(l — 2M/r) + const is 
the Regge- Wheeler tortoise coordinate while the poten- 
tial V{r) is given by 



V(r) 



-2M 



£(£ + 1) 2M 



(3) 



(ii) as any physical wave, must have a purely ingoing 
behavior at the event horizon at r = 2M, i.e. 



®l(r) ~ e" 

r* — *— oo 



(4) 



(iii) At spatial infinity r 
behavior 



-oo, has the asymptotic 
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-iwr*+i^7r/2 



S({uj) 



-\-itur* —Htt/2 



(5) 

For certain complex values of to, both Se and Tg have 
a simple pole but Se/Tg is regular. These values are 
the frequencies of the QNMs, which we can define (see 
Eqs. J3J and iJSJ) as the solutions of the wave equation 
which represent a purely outgoing wave at infinity and 
a purely ingoing wave at the horizon. We now denote 
by LUi„ = <Jf^ — iTe n /2 with Jf^ > and Tg n > the 

QNM frequencies, u>^ representing the frequency of the 
oscillation corresponding to the QNM and Tg n represent- 
ing its damping. In the immediate neighborhood of coi n , 
Se(oj) has the Breit-Wigner form, i.e., 



Se{uj) oc 



T ep /2 



(6) 



Using the CAM method, we can provide a physical 
picture of the scattering process in term of diffraction by 
surface waves and a physical explanation of the mecha- 
nism of QNM excitation valid for high frequencies. By 
means of a Watson transformation 1] applied to the scat- 
tering amplitude Q, we can write Q 



A(5 A _ 1/2 (w)-l) 



cos7rA 



P\-i/2(-cos0) dX. 

(7) 
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FIG. 1: The Watson integration contour. 



Here C is the integration contour in the complex A-plane 
[lj illustrated in Fig. 2] The Watson transformation 
permits us to replace the ordinary angular momentum 
£ by the complex angular momentum (CAM) A. Here 
S\-i/2( LU ) is now an analytic extension of Se(u>) into the 
complex A-plane which is regular in the vicinity of the 
positive real A axis. Moreover, Pa-i/2( 2 ) is the hyperge- 
ometric function F(-X + 1/2, A + 1/2; 1; (1 - z)/2). 

We can then deform the path of integration in Eq. J7J| , 
taking into account the possible singularities. The only 
singularities that are encountered are the poles of the 
S-matrix lying in the first quadrant of the CAM plane 
• They are known as Regge poles 0, Q an d we shall 
denote them by X n (u>), the index n — 1,2,... permit- 
ting us to distinguish between the different poles. By 
Cauchy's Theorem we can then extract from Eq. Q the 
contribution of a residue series over Regge poles given by 



-in ^ X n (uj)r n (uj) 



LU ^— ' cos (nXn(ui)) 



P 



(— COS^ 



(8) 



where r n (cu) = residue (<Sa-i/2(w)) a _ a , It should be 
noted that / differs from fp by a background integral 
P| which does not play any role in the resonance phe- 
nomenon. By using the asymptotic expansion 



Pa-i/2(-cos< 



e iA(7r-0)-i7r/4 _|_ e -iA(7r-0)+i7r/4 

(27rAsin6') 1 /2 



as |A| — > oo, valid for |A| sm9 > 1, as well as 



cos7rA " ' 
m=0 



+ OC 

2i e " r(2m+1)(A ~ 1/2) 
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which is true if Im A > 0, we can write 
2tt \ n {uj)r n (u) 



n = 

o 

m=0 



[2^A„Hsinfl] 1/2 

+ 00 



iX n (lj) (0+2m7r) — im 77+271-/4 



*A Tl (aj)(27r — 0+2m7r)— zm7r— 27r/4 



.(9) 



In Eq. (2J, terms like exp[iA n (w)(#)] and exp[iA„(o;)(27r — 
0)] correspond to surface wave contributions. Because 
a given Regge pole X n (uj) lies in the first quadrant of 
the CAM plane, exp[iA„(w)(#)] (resp. exp[iA„(u;)(27r — 
9)]) corresponds to a surface wave propagating counter- 
clockwise (resp. clockwise) around the black hole and 
Re X n (u>) represents its azimuthal propagation constant 
while Im X n (u>) is its damping constant. The exponential 
decay exp[— Im X n (uj)8] (resp. exp[— Im A„(w)(27r — 9)]) 
is due to continual rcradiation of energy. Moreover, in 
Eq. @, the sum over m takes into account the multiple 
circumnavigations of the surface waves around the black 
hole as well as the associated radiation damping. Finally, 
the presence of the factor exp[— imir] in Eq. © should 
be also noted: it accounts for the phase advance due to 
the two caustics on the scattering axis. 

The resonant behavior of the black hole can now be 
understood in terms of surface waves. As u> varies, each 
Regge pole X n (uj) describes a Regge trajectory Q in the 
CAM plane. When the quantity Re X n (oj) coincides with 
a half-integer (a half-integer but not an integer because of 
the caustics), a resonance occurs. Indeed, it is produced 
by a constructive interference between the different com- 
ponents of the n-th surface wave, each component cor- 
responding to a different number of circumnavigations. 
Resonance wave frequencies are therefore obtained 
from the Bohr-Sommerfeld type quantization condition 

'"'i f+\ £ = 0,1,2,.... (10) 

By assuming that u> is in the neighborhood of and 
using Re A„(a->) ^ Im X n (cj) (which can be numerically 
verified, except for very low frequencies), we can expand 

X n (u>) in a Taylor series about , and obtain 



Re A, 



X n (u) W £+ - + 



d ReA„(w) 



-i Im A n (w^° ) ). 



(ii) 

Then, by replacing Eq. (|11|) in the term cos (7rA T j(o;)) 
of Eq. (JSJ), we show that fp(w,9) presents a resonant 
behavior given by the Breit-Wigner formula JHJ with 



Im X n (u) 



d Re X n (u)/du> 



(12) 



Eqs. I|10[l and l|12f) are kind of semiclassical formulas 
which permit us to determine the location of the reso- 
nances from Regge trajectories. 
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FIG. 2: The Regge poles A n (w) followed for u = 
n — 1, 2, 3, 4 for a scalar perturbation (2M = 1). 
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We now discuss the numerical aspects of our work. In 
order to determine the location of the QNM frequencies 
from Eqs. (|10fl and (|12(l . we need the Regge trajectories. 
For a given lu real and positive, we search for a complex 
value of I such that the coefficient 1/Ti(u>) is zero but 
Si(lo) /Ti(lv) is not. Such a value is then a pole of Se(u>) 
and provides immediately the value of the associated pole 
for S\-i/2{u))- Leaver [12j presents a method for finding 
the zeros of 1/T^(w). Though he uses his method to 
find the QNM frequencies, his method is equally valid, 
mutatis mutandis, for finding the Regge poles {A n = £ n + 
1/2} of Si(oj). Leaver writes the solution to Eq. as 
an infinite sum 



fc=0 



r-2M 



(13) 



where p = —i2Mu>. The recurrence relations between the 
ak are given by Leaver: 



where 

ak = 
(3 k = 

7fc = 



a k a k +i + f3 k a k + 7fc fl fc-i = 0, 



k 2 + 2fc(p+l) + 2p+l, 



(14) 



(15) 

-[2k 2 + 2k(4p + 1) + 8p 2 + 4p + l{l + 1) + 1], 

(16) 

2 + 4kp + 4p 2 . (17) 



Leaver noted that the coefficient 1/Tt(u>) has a zero 
whenever the sum ^ ak converges. He translates this 
requirement into an infinite-fraction equation involving 
the coefficients a, (3 and 7. Majumdar and Panchapake- 
san gave an alternative condition using the Hill determi- 
nant method ^4| : a set of parameters {I, lj} that give a 
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TABLE I: A sample of QNM frequencies for scalar perturba- 
tions (2M = 1). 







Exact" 


Exact" 


Semiclassical 


Semiclassical 


t 


n 


/ >(°) 


r^n/2 


, ,(°~> 


Efn/2 





1 


0.22091 


-0.20979 


0.19009 


-0.23120 




2 


0.17223 


-0.69610 


0.09743 


-0.33189 




3 


0.15148 


-1.20216 


0.09695 


-0.31452 


1 


1 


0.58587 


-0.19532 


0.58336 


-0.19653 




2 


0.52890 


-0.61252 


0.51960 


-0.59418 




3 


0.45908 


-1.08027 


0.44000 


-0.93961 


2 


1 


0.96729 


-0.19352 


0.96669 


-0.19371 




2 


0.92716 


-0.59121 


0.92504 


-0.58871 




3 


0.86109 


-1.01712 


0.85890 


-0.97714 



"from Ref.lll 



convergent sum solve 



D = 



7i 



a 

72 



a 2 



= 0. 



(18) 



This technique has the advantage of being more numeri- 
cally tractable than Leaver's, but it might be only valid 
for the Regge poles that lie close to the real axis of the 
CAM plane. Both Leaver's original method and the Hill 
determinant method have been applied to the search for 
Regge poles in our study. We found excellent agree- 
ment between Leaver's method and the Hill-determinant 
method for small values of Regge mode index n. We 
also agreed with Andersson's values given in Ref. @- 
Fig. 12 exhibits the Regge trajectories numerically calcu- 
lated while Table[I]presents a sample of QNM frequencies 
calculated from the semiclassical formulas (|10ll and (|12|l . 
A comparison between the "exact" and the semiclassical 
spectra shows a good agreement, except for very low fre- 
quencies. Furthermore, the semiclassical theory permits 
us to classify the resonances in distinct families, each 
family being associated with one Regge pole and there- 
fore to understand the meaning of the indices n and I 
introduced to denote the QNM frequencies iOi n . 

At large u, the position of the Regge poles very closely 
adheres to the asymptotic form: 

A„(w) w 3V3Mu + i(n- 1/2) n = 1,2,3,... (19) 

Eq. i|19|) leads us to conclude that the n-th surface wave 
is localized near the unstable circular photon orbit at 
R = 3M by the following consistency argument: by rein- 
stating dependence on Schwarzschild time t into Eq. 10, 
we can see that the surface waves circle the black hole in 
time T — 2nKe X n (ui)/u> 2ir3\^3M for large u>. Fur- 
thermore, a photon on the circular orbit with constant 
radius R takes the time T' = 2nR/(l - 2M/R) 1/2 to cir- 
cle the black hole (this result can be found by integrating 
the Schwarzschild metric ds 2 = 0). By equating T and 
T', we obtain R = 3M. Moreover, by using Eqs. (|10fl . 



(|12fl and l|19|) , we recover the well-known high-frequency 
behaviors (see Refs. and0 and references therein) 



(o) 
In 



1/2 



3\/3M 



and 



1/2 



3V3M 



(20) 



These asymptotic behaviors have been obtained under 
the hypothesis Re A„(w) » Im X n (u)) and are therefore 
valid for I ^> n. 

Our approach also applies to electromagnetic and grav- 
itational scattering. In Fig. [3] and Table [H] we present 
some results for the latter case. A comparison between 
the exact and the semiclassical spectra shows a rather 
good agreement. 

Finally, it seems to us necessary to emphasize that the 
CAM method we have developed here is based on two 
assumptions about the Regge poles A„(oj): They must 
formally satisfy |A„(w)| — > +oo as well as Re X n (uj) ^> 
Im A„(w) as u — > oo. As a consequence, our method is a 
"high-frequency" approach which permits us to recover 
and to physically interpret the spectrum of the QNM 
frequencies lying in the region 2M Re u> > 0.2 of the 
complex w-plane. It should be noted that the highly 
damped QNM (for gravitational perturbations) as well 
as their associated frequencies which are connected with 
the area spectrum of the black hole (see for the first 
paper on this subject and 0] and references therein for 
recent works in this domain) can neither be understood in 
the framework of our approach nor interpreted in terms 
of surface waves orbiting around the black hole near the 
unstable circular orbit at r = 3M: Indeed, when the 
index n is very large, the frequencies for these highly 

damped QNM are asymptotically given by momna 



,(°) 



In 3 
8ttM 



and 



1/2 



AM 



(21) 



and they therefore lie in the region 2M Re u> <C 0.2 and 
2M |Im u\ ^> 1 of the complex w-plane. 



III. CONCLUSION 

To conclude, we would like first to comment on some 
aspects of our work and then to consider some possible 
extensions of the CAM approach in gravitational physics: 

- We have established the connection between Anders- 
son's surface waves and the QNM complex frequencies 
of the Schwarzschild black hole in the particular con- 
text of wave scattering. This connection is more gen- 
eral and could be also obtained directly from the Regge- 
Wheeler equation, by extending the CAM method devel- 
oped by Sommcrfeld |2fJ as an alternative to the Watson 
approach 

- Because the gravitational radiation created in 
many black hole processes is dominated at intermedi- 
ate timescales by QNMs, it can be always interpreted 
in terms of surface waves. Whatever the perturbation 
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2 4 6 8 10 12 14 

FIG. 3: The Regge poles A n (tj) followed for ui = — > 5, 
n = 1,2, 3, 4, for gravitational perturbations (2M = 1). 



cently proposed in the context of superstring theory |2l| . 
With this aim in view, the (2+l)-dimensional BTZ black 
hole seems to us very interesting because, in that 
particular space-time, the wave equation can be solved 
exactly (see, e.g., |23|) and therefore Regge poles can be 
analytically obtained. 

(iv) The study of artificial black holes |24| . Here, it would 
be possible to benefit from the formidable CAM machin- 
ery developed in electromagnetism, optics and acoustics. 
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TABLE II: A sample of QNM frequencies for gravitational 
perturbations (2M = 1). 







Exact" 


Exact" 


Semiclassical 


Semiclassical 


£ 


n 




Efn/2 




Tfn/2 


2 


1 


0.74734 


-0.17793 


0.75812 


-0.17644 




2 


0.69342 


-0.54783 


0.78134 


-0.49976 




3 


0.60211 


-0.95655 


0.77683 


-0.83126 


3 


1 


1.19889 


-0.18541 


1.20332 


-0.18455 




2 


1.16529 


-0.56260 


1.20089 


-0.54409 




3 


1.10337 


-0.95819 


1.18284 


-0.89846 


4 


1 


1.61836 


-0.18832 


1.62052 


-0.18792 




2 


1.59326 


-0.56866 


1.61119 


-0.56020 




3 


1.54542 


-0.95982 


1.58849 


-0.92907 



"from Chapter 4 of Ref.0 or Refs.0 andU 



which modifies the geometry of a Schwarzschild black 
hole, it leads to the excitation of surface waves localized 
close to the unstable photon orbit. During its repeated 
circumnavigations, a given surface wave decays by ra- 
diating away its energy, leading to the damped ringing 
of the geometry. This mechanism is analogous to the 
damped ringing of the Earth which occurs during several 
days after a large earthquake and which can be explained 
in term of the Rayleigh surface wave. 

- The CAM method could be naturally introduced in 
many other areas of the physics of relativistic stars and 
black holes. Such a program could include in particular 
the following topics: 

(i) The interpretation of the QNM of neutron stars and 
of the Kerr black hole in terms of surface waves. For the 
Kerr black hole, such an approach could easily explain the 
splitting of the quasinormal frequencies due to rotation. 

(ii) The analysis of Hawking radiation as well as of Kerr 
black hole superradiance. 

(iii) The study of black holes immersed in asymptotically 
anti-de Sitter space-times. The CAM approach could 
provide new tests of the AdS/CFT correspondence re- 
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